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ABSTRACT 

This paper presents the results from the first wide and deep dual narrow-band survey 
to select Ha and [Oil] line emitters at z = 1.47 ± 0.02, exploiting synergies between the 
United Kingdom InfraRed Telescope and the Subaru telescope by using matched narrow- 
band filters in the H and z' bands. The Ha survey at z = 1.47 reaches a 3 a flux limit of 
Fhc* ~ 7 x 10 _17 ergs _1 cm -2 (corresponding to a limiting SFR in Ha of « 7M yr _1 ) 
and detects « 200 Ha emitters over 0.7 deg 2 , while the much deeper [Oil] survey reaches an 
effective flux of « 7 x 10" 18 erg s _1 cm -2 (SFR in [On] of ~ 1 M yr" 1 ), detecting « 1400 
z = 1.47 [Oil] emitters in a matched co-moving volume of ~ 2.5 x 10 5 Mpc 3 . The combined 
survey results in the identification of 190 simultaneous Ha and [Oil] emitters at z = 1.47. 
Ha and [Oil] luminosity functions are derived and both are shown to evolve significantly 
from z ~ in a consistent way. The star formation rate density of the Universe at z = 1.47 is 
evaluated, with the Ha analysis yielding psfr = 0.16 ± 0.05 M yr _1 Mpc -3 and the [Oil] 
analysis Psfr = 0.17 ± 0.04 M yr _1 Mpc -3 . The measurements are combined with other 
studies, providing a self-consistent measurement of the star formation history of the Universe 
over the last ~ 11 Gyrs. By using a large comparison sample at z ~ 0.1, derived from the 
Sloan Digital Sky Survey, [Oll]/Ha line ratios are calibrated as probes of dust-extinction. 
Ha emitters at z ~ 1.47 show on average Ah« ~ lmag, the same as found by SDSS in 
the local Universe. It is shown that although dust extinction correlates with SFR, the relation 
evolves by about ~ 0.5 mag from z ~ 1.5 to z ~ 0, with local relations over-predicting the 
dust extinction corrections at high-z by that amount. Stellar mass is found to be a much more 
fundamental extinction predictor, with the same relation between mass and dust-extinction 
being valid at both z ~ and z ~ 1.5, at least for low and moderate stellar masses. The 
evolution in the extinction- SFR relation is therefore interpreted as being due to the evolution 
in median specific SFRs over cosmic time. Dust extinction corrections as a function of optical 
colours are also derived and shown to be broadly valid at both z ~ and z ~ 1.5, offering 
simpler mechanisms for estimating extinction in moderately star-forming systems over the 
last ~ 9 Gyrs. 

Key words: galaxies: high-redshift, galaxies: luminosity function, cosmology: observations, 
galaxies: evolution. 



* Based on observations obtained using both the Wide Field CAmera on 
the 3.8m United Kingdom Infrared Telescope (UKIRT), as part of the High- 
z Emission Line Survey (HiZELS), and Suprime-Cam on the 8.2-m Subaru 
Telescope, which is operated by the National Astronomical Observatory of 
Japan. 



1 INTRODUCTION 

It is now clear that the "epoch" of galaxy formation occurs at z > 1, 
as surveys measuring the star formation rate density (psfr) as a 
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function of epoch show that psfr rises steeply out to z ~ 1 (e.g. 
|Lilly et al.|[T^96l |Hopkins & Beacom|2006T >, but determining the 
redshift where psfr peaked at z > 1 is more difficult. 

Accurately determining psfr requires the selection of clean 
and well-defined samples of star-forming galaxies over representa- 
tive volumes. In practice, such selection is done by detecting sig- 
natures of massive stars (being very short-lived, their presence im- 
plies recent episodes of star formation). The high luminosities of 
such massive stars allow them to be traced up to very high red- 
shift, and to estimate star formation rates (SFRs). Such stars emit 
strongly in the ultra-violet (UV), and although this is often sig- 
nificantly absorbed, it is then re-emitted through a variety of pro- 
cesses, resulting in detectable signatures such as strong emission 
lines or thermal infra-red emission from heated dust (see discus- 
sion in Hopkins et al. 2003). Ideally, the use of different tracers of 
recent star formation would provide a consistent view, but in reality, 
because they have different selection biases, and require different 
assumptions/extrapolations, measurements are often significantly 
discrepant (e.g. Hopkins & Beacom 2006). A major issue is the 
difficulty and uncertainty in correcting for (dust) extinction, espe- 
cially for UV and optical wavelengths, leading to large systematic 
uncertainties (and potentially missing entire populations). An addi- 
tional complication comes from surveys at different epochs making 
use of different indicators due to instrumentation and detection lim- 
itations. For example, while for < z £ 0.4 the evolution of psfr 
is typically estimated using Ha luminosity, because it is easily tar- 
geted in the optical (e.g. |Ly et al.|2007| |Shioya et al.|2008[ [Dale] 
|et al.|2010] >, at higher redshifts the line is redshifted into the near- 
infrared, and [Oil] 3727 - hereafter [Oil] - luminosity (much more 
affected by extinction than Ha, and also a metallicity-dependent 
indicator) is used instead (e.g. Zhu et al. 2009 ; Bayliss et al. 201 1). 

The [Oil] emission line can be traced in the optical window up 
to z ~ 1.5 and many authors have attempted to measure the evolu- 
tion of [Oil] luminosity function up to such look-back times (e.g. 



Teplitz et al.||2003| 


Hogg et al.|1998||Gallego et al.||2002| |Taka-| 


|hashi et al.|2007||Z1 


m et al.|2009|>. Such studies have identified a 



|Zhu et al.|2009| ), although the evolution seems to be slightly dif- 
ferent from that seen in the Ha luminosity function. Part of the 
differences may well arise from the difficulty in using [Oil] lumi- 
nosity density directly as a star-formation rate density indicator: 
[Oil] luminosity is calibrated locally using Ha (see e.g. [Aragon- 
Salamanca et al. 2003, Mouhcine et al. 2005 for a comparison 
between both emission-lines at low-z), but the actual calibration 
depends on metallicity and dust-extinction (c.f. |Jansen et al.j 2001 ; 
Kewley et al. 2004). These issues have been reported by many stud- 
ies, mostly using data from the local Universe or low redshift (see 
e.g. |Gilbank et al. 2010 for a comparison between both lines and 
other indicators in the Sloan Digital Sky Survey), for which it is 
possible to measure both Ha and [Oil]. However, so far it has not 
been possible to directly compare the Ha and [Oil] indicators at 
z > 1 using large-enough samples to test, extend and improve our 
understanding. 

While some Ha studies beyond z ~ 0.4 have been conducted 
in the 90's (e.g. Bunker et al.|1995[|Malkan et al.|1995| >, tracing the 
Ha emission line in the infrared, the small field-of-view of infrared 
detectors at that time made it very hard to detect more than ^1 — 2 
emitters at z > 1 (e.g. jvan der Werf et al. 2000) using ground-based 
narrow-band surveys. Slitless spectroscopy using NICMOS on the 
Hubble Space Telescope (HST) provided a space-based alternative 
to make significant progress, particularly by allowing much deeper 
observations than those from the ground (e.g. McCarthy et al. 1999 ; 



|Yanetal.| 1999, Hopkins et al. 2000). Further progress is now be- 
ing obtained with WFC3 on HST, as recent studies (e.g. Ate k et ah] 
|2010[|Straughn et al.|20lT]|van Dokkum et al.|201 1\ take advantage 
of the increased sensibility and the wider field-of-view (although 
still relatively small when compared to wide-field ground-based in- 
frared detectors) of WFC3 to simultaneously look for fainter emit- 
ters and increase the sample sizes. Furthermore, the development 
of wide-field near-infrared cameras in the last decade has finally 
made it possible to conduct ground-based large area Ha surveys 
(necessary to overcome cosmic variance) which can look for such 
emission-line galaxies all the way up to z ~ 2.5. 

HiZELS, the High-redshift(Z) Emission Line Surve^ (Geach 



|eFaLl|2008] |Sobral et aL]|2009a| hereafter S09) is a Campaign 
Project using the Wide Field CAMera (WFCAM) on the United 
Kingdom Infra-Red Telescope (UKIRT) and exploits specially- 
designed narrow-band filters in the J and H bands (NB j and NBh), 
along with the H2SI filter in the K band, to undertake panoramic, 
moderate depth surveys for line emitters. HiZELS is primarily tar- 
geting the Ha emission line redshifted into the near-infrared at 
z = 0.84, z = 1.47 and z = 2.23. The main HiZELS survey aims 
to cover « 7 deg 2 (to overcome cosmic variance) and is detecting 
~ 1000 emitters over volumes of > 10 6 Mpc 3 with each filter, 
reaching limiting Ha (observed) SFRs of « 3, 6 and 30 M yr _1 
at z = 0.84, z = 1.47 and z = 2.23. Such data will pin down 
the likely peak of /?sfr and provide detailed information about the 
population of star-forming galaxies at each epoch (see Best et al. 
|20T0] ). 

HiZELS provides precise measurements of the evolution of 
the Ha luminosity function from z = 0.0 to z = 2.23 (Geach et al. 
2008; S09), while the contribution from much deeper Ha surveys 
over smaller areas (e.g. |Hayes et"aL]|2010| |Ly et al.||2011| ) offers 
the possibility of also exploring the faint-population. The Ha lumi- 
nosity function evolves significantly, mostly due to an increase by 
more than one order of magnitude in L^ a (S09), the characteristic 
Ha luminosity, from the local Universe to z = 2.23. In addition, 
|Sobral et ah] ( |2011| > found that at z = 0.84 the faint-end slope of 
the luminosity function is strongly dependent on the environment, 
with the Ha luminosity function being much steeper in low den- 
sity regions and much shallower in the group/cluster environments. 
Whether this is also found at higher redshifts remains unknown. 

Indeed, while it is now possible to detect Ha emission over 
very wide areas up to 2 > 2, distinguishing between Ha and 
any other emission line at any redshift is a challenge, particu- 
larly over large areas where it is unfeasible to assemble ultra- 
deep multi- wavelength data over dozens of bands. Matched dual 
narrow-band surveys offer a solution to the problem; in particu- 
lar, since the NB921 narrow-band filter on Subaru is able to probe 
the [Oil] emission line for the same redshift range (z ~ 1.47) as 
the HiZELS narrow-band H filter on UKIRT probes Ha (see Fig- 
ure [TJ, a matched and sufficiently deep survey would not only pro- 
vide a simple, clean selection (z = 1.47 sources will be detected 
as emitters in both data- sets), but also provide a means of directly 
comparing Ha and [Oil] at z ~ 1.5 for large samples for the first 
time. 

This paper presents deep narrow-band imaging using the NBh 
filter at A = 1.617|xm, as part of HiZELS, over 0.79 deg 2 in the 
UKIDSS Ultra Deep Survey (Lawrence et al.||2007] ) field (UDS) 
and combines the data with ultra-deep NB921 imaging taken us- 



1 For more details on the survey, progress and data releases, see 
http://www.roe.ac.uk/ifa/HiZELS/ 
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Table 1. Observation log for the NB H observations of the UDS field, taken using WFCAM on UKIRT, during 2008, 2009 and 2010. Limiting magnitudes 
(3 a) are the average of the four frames which constitute a field, based on the the measurements of 10 6 randomly placed 2" apertures in each frame. Integration 
times given in brackets present the total integration times obtained at the telescope prior to rejection of data taken in poorer conditions (see §2.1} . 
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Figure 1. A comparison between the NBh redshift coverage of the Hen 
emission line and the NB921 redshift coverage of the [On] emission-line, 
revealing that the two filters are extremely well matched and ideal for a dual 
Hen- [On] narrow-band survey. The [On] coverage is only slightly wider (in 
redshift) than the Ha coverage - but this is shown to have little to no effect 
on the analysis (c.f. Section [43) . 

ing Suprime-Cam on the Subaru telescope, to explore an extremely 
well-matched Ha- [Oil] narrow-band survey over 0.70 deg 2 . 

The paper is organised as follows: §2 presents the obser- 
vations, data reduction, source extraction and catalogue produc- 
tion. In §3, emission-line galaxies are selected and, using colours 
and photometric redshifts, the samples of Ha and [Oil] selected 
z = 1.47 emitters are presented. §4 presents the Ha and [Oil] lu- 
minosity functions and the derived star-formation rate density at 
z = 1.47, together with an accurate measurement of their evolu- 
tion. In §5, the large sample of robust Ha- [Oil] emitters is used to 
conduct line ratio studies and compare these with a large sample 
selected locally from the Sloan Digital Sky Survey (SDSS). These 
are used to investigate the extinction properties of the galaxies, and 
how these evolve with redshift. Finally, §6 outlines the main con- 
clusions. An H = 70kms _1 Mpc -1 , VLm — 0.3 and = 0.7 
cosmology is used. Magnitudes are in the Vega system, except if 
noted otherwise. 



2 DATA AND SAMPLES 

2.1 Near-infrared NB H imaging with UKIRT 

The UKIDSS UDS field was observed with WFCAM on UKIRT 
using a set of custom narrow-band H filters (NBh, A = 1.617 \xm, 
AX = 0.021 um), as detailed in Table 1. WFCAM's standard 



"paw-print" configuration of four 2048 x 2048 (0.4" pixel -1 ) de- 
tectors offset by ~ 20' was mac rostepped four time s to cover a con- 
tiguous region of ~ 55' x 55' ( [Casali et al.|2007) , with individual 
narrow-band exposures of 100 s. The conditions were mostly pho- 
tometric. The Non Destructive Read (NDR) mode was used for all 
narrow-band observations to minimise the effects of cosmic rays in 
long exposures. The observations were obtained by jittering around 
14 different positions in each of the 4 pointings. The UKIDSS UDS 
if -band image (24 AB, 5a) overlaps entirely with the full narrow- 
band image, yielding a total survey area of 0.79 deg 2 . 

A dedicated pipeline has been developed for HiZELS 
(PfHiZELS, c.f. S09 for more details), and was used to fully re- 
duce the data. Very briefly, the pipeline dark- subtracts the data and 
median combines dark- subtracted data (without offsetting) to ob- 
tain a first-pass flat field for each jitter pattern. The latter is used to 
produce a badpixel mask for each chip by flagging pixels which de- 
viate by more than 5cr from the median value. Frames are individ- 
ually flattened and individual source masks are produced for each 
frame. A new final flat is then created masking detected sources 
(for each jitter pattern) and all data are flattened. A world coordi- 
nate system is fitted to each frame by querying the USNO A2.0 
catalogue (typically ^ 70 — 80 stars per frame). Finally, all indi- 
vidual reduced frames were visually inspected, resulting in the re- 
jection (for each of the four cameras) of 12 frames for the NW field 
and 8 frames in the SW field due to bad quality (FWHM> 1.1" 
or Ah > 0.25 mag). Final frames which passed the data-control 
were de-jittered and median co-added with SWarp (Bertin 1998), 
including a background mesh-based sky subtraction optimized for 
narrow-band data, and masking of bad-pixels and cosmic-rays. It 
should be noted that WFCAM frames suffer from significant cross- 
talk artifacts (toroidal features at regular pixel intervals from bright 
sources in the read-out direction); because these are linked to a 
"physical" location, they can only be removed from the source cat- 
alogue (see j ]2.2| l. 

Narrow-band images were photometrically calibrated (inde- 
pendently) by matching ~ 70 stars with ttih = 11 — 16 per frame 
from the 2MASS All-Sky catalogue of Point Sources ( |Cutri et aL] 
|2003| ) which are unsaturated in the narrow-band frames. 

2.2 NBh Source Extraction and Survey Limits 

Sourc es were extracted using S EXTRACTOR ( [Bertin & Arnouts| 
1996 ). Photometry was measured in apertures of 2" diameter. In 
order to clean spurious sources from the catalogue (essential to re- 
move cross-talk artifacts), the final images were visually inspected; 
this revealed that sources brighter than ~ 11.5 mag (NBh, Vega) 
were surrounded by a number of artifacts detected within "bright 
halos", as well as cross-talk. Fainter sources (up to 16 mag) showed 
only cross-talk features. Sources fainter than 16 mag did not pro- 
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duce any detectable cross-talk at the depth of the observations. 
Cross-talk sources and detections in the halo regions were removed 
from the catalogue separately for each frame, which greatly simpli- 
fies their identification. 

The average 3 a depth of the entire set of NBh frames is 
21.2 mag; this is measured using a set of 10 6 randomly placed 
apertures per frame. Above the 3 a threshold in each frame, the 
narrow-band imaging detects a total of 23394 sources (5904, 4533, 
6946 and 601 1 in the NE, NW, SE and SW pointings, respectively) 
across 0.78 deg 2 (after removal of regions in which cross-talk and 
other artifacts caused by bright objects are located) 



2.3 Optical NB921 imaging with Subaru 

Archival Subaru/Suprime-Cam NB921 data of the UDS field are 
available, taken by |Ouchi et al.| ( [20091 |20T0] l. The field was ob- 
served with Suprime-Cam on Subaru as part of the Subaru/XMM- 
Newton Deep Survey (SXDS, |Ouchi et al.1[2008] > during 2005- 
2007. The NB921 filter is centered at 9196 A^ith a FWHM of 
132 A, and was used to cover the field with 5 pointings with a to- 
tal integration of 45.1 hours (individual exposures ranging from 
8 to 10 hours) - see |Ouchi efaLl ( [20T0] >. The raw NB921 data 
were downloaded from the archive and reduced with the Suprime- 
Cam Deep field REDuction package (SDFRED, [Yagi et al.|2002l 
|Ouch i et al. 2004) and IRAF. The combined images were aligned to 
the public z'-band im ages of Subaru-XMM Deep Survey (SXDS; 
|Furusawa et"aLl|2008| ) and PSF matched (FWHM= 0.8"). The 
NB921 zero points were determined using z data, so that the z - 
NB921 colour distributions of SXDS would be consistent with that 
of Subaru Deep Field data (Kashikawa et al. 2004). Source detec- 
tion and photometry were performed using SEXTRACTOR ([Bertin 
& Arnouts 1996). Sources were detected on each individual NB921 
image and magnitudes measured with 2" diameter apertures. The 
average NB921 3a limiting magnitude is estimated to be 26.3 (AB) 
by randomly placing 10 6 2" apertures in each frame; down to that 
depth, 347341 sources are detected over the entire 1 deg 2 area (5 
Suprime-Cam pointings). Note that NB921 and z' magnitudes are 
given in the AB system. 



3 SELECTION 

3.1 Narrowband excess selection 

Potential line emitters (NBh and NB921) are selected according to 
the significance of their broad-band — narrow-band (BB — NB) 
colour, as they will have (BB — NB) > 0. However, true emitters 
need to be distinguished from those with positive colours due to 
the scatter in the magnitude measurements; this is done by quanti- 
fying the significance of the narrowband excess. The parameter E 
(see S09) quantifies the excess compared to the random scatter ex- 
pected for a source with zero colour, as a function of narrow-band 
magnitude (c.f. S09). 

Neither of the narrow-band filters falls at the center of the 
broad-band filters and thus objects with redder colours will tend 
to have a negative (H— NBh) colour, and a positive (z ; -NB921) 
colour, while bluer sources will have (if— NBh) > and (z - 
NB921) < 0. This will not only affect the selection of emission line 
objects, but can also result in an over/under-estimation of emission- 
line fluxes because it will lead to an under/over-estimation of the 
continuum. Fortunately, this can be broadly corrected by consider- 
ing the broad-band colours of each source, as in S09. For the NBh 



H-NBh 



NB921 




Figure 3. The top panels show a strong NB H emitter, clearly revealed after 
subtracting the H continuum. The bottom panels show the same sky area 
(~ 30" x 30") in NB921 and z' , revealing that the NB H emitter is also 
a strong NB921 emitter (and thus a Hq;+[Oii] z ~ 1.47 emitter). Note 
the depth difference, and how non-emitters disappear when the continuum 
estimated flux is subtracted from the narrow-band flux (other fainter NB921 
emitters are also revealed). 



data, this is done by studying (H— NBh) as a function of (J — H) 
colour for all sources detected in the NBh frames. A linear fit is 
derivecj^jand is used to correct the initial H magnitudes to produce 
an effective H' magnitude appropriate for estimation of the contin- 
uum contribution at the wavelength of the NBh filter. This assures 
a mean zero (H' '— NBh) as a function of (J — H), and also results 
in no significant trend as a function of (H — K). The correction is 
given by: 



(H' - NB t 



(1) 



(H - NB H ) - 0.1(J -H) - 0.03. 

For the NB921 data, a similar approach is taken, but using z — J 
colours insteacf^] Colours are corrected by using the best linear fit: 

(z" - NB921) = (z - NB921) - 0.05(z' - f) + 0.15. (2) 



Emission line fluxes, Fii ne , 
puted using: 

Inb 



and equivalent widths, EWii ne , are com- 



Fline = AAj> 



1 - (AX N b/AXbb) 



and 



EWii 



AAj\ 



Inb — Sbb 



(3) 



(4) 



fBB — /nB^Xnb/AXbb) 

where AXnb and AXbb are the FWHMs of the narrow- and 
broad-band filters (21 1 A and 132 A for NB H and NB921 ; 2893 A 
and 955 A for H and z respectively), and /nb and fsB are the 
flux densities measured for the narrow and broad-bands, respec- 
tively. Note that fs b is computed by using the corrected BB' mag- 
nitudes, which are a much better approximation of the continuum 
for the present purposes and guarantee that the median flux dis- 
tribution is zero. Note that the latter formula is only valid because 
magnitudes and colours have been corrected to guarantee a flux and 
colour distribution centered at zero. 

The selection of emission-line candidates is done following 

2 Only sources within ±2<j of the general scatter around the median 
H— NBh are used. Furthermore, in order to improve the fit for galaxies, 
potential stars (selected as sources satisfying B — z > 13.5(H — K) + 2.0, 
AB) are excluded and a separate fit is done just for the latter. 

3 For sources with no clear J detection, a statistical correction of +0.05 is 
applied. 
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NB H (mag) NB921 (mag, AB) 

Figure 2. Narrow-band excess as a function of narrow-band magnitude for NBh (left panel, Vega magnitudes) and NB921 (right panel, AB magnitudes) data. 
These show > 3a detections in narrow-band imaging and the lines present the average 2.5 E and 3.0 E colour significances for NBh and NB921, respectively 
(for the average depth, but note that the analysis is done individually for each frame). Also, due to the very high number of detected sources in NB921, only 
1 in every 5 sources is shown. The horizontal dashed lines present the equivalent width cuts used for NBh and NB921 data - these correspond to z = 1.47 
rest-frame EW limits of 20 A for Hen and 10 A for [On]. Narrow-band excess sources in both NBh and NB921 (double-emitters) are shown, and constitute a 
robust sample of Hcn/[Oii] candidate emitters - see Figure[3]for a visualisation of one such emitters. 



S09 (see also |Ouchi eFaTHIOTOl ). Narrow-band sources in the H 
band (NBh) are selected as line emitters if they have an individual 
> 3 a detection in NBh and present a colour excess significance 
of E > 2.5 (which broadly corresponds to a flux limit, see S09), 
and an observed EW > 50 A (corresponding to EW > 20 A rest- 
frame for the Ha line at z = 1.47). NB921 sources are selected 
as potential emitters if they have an individual > 3 a detection in 
NB921 and present a colour excess significance of E > 3.0 (as 
the data are much deeper, a higher significance cut can be applied) 
and EW> 25 A (corresponding to EW > 10 A for [Oil] emitters 
at z = 1.47). The EW cuts are applied to avoid including bright 
foreground objects with a large significance and a steep contin- 
uum across the H or z bands, and were chosen to reflect the gen- 
eral scatter around the zero colour at bright magnitudes (thus the 
difference between NBh and NB921) and to allow a good selec- 
tion of both Ha and [Oil] emitters. Figure [2]presents the corrected 
broad-band — narrow-band colours as a function of narrow-band 
magnitude, including the selection criteria and the sample of NBh 
and NB921 emitters. Figure [3]presents examples of emitters drawn 
from the samples (for each filter). 



3.2 The samples of narrow-band emitters 

Narrow-band detections below the estimated 3 a threshold were not 
considered. Due to the depth of the UKIDSS-UDS data (DR5), all 
extracted NBh sources have > 3 a broad-band H detections (down 
to ~ 24.9 AB); for the much deeper NB921 data, only sources with 
> 3a detections in z' (z f < 26.7) are considerecjj The average 
3 a line flux limit is 7 x 10~ 17 erg s _1 cm -2 for the NBh data 



4 This results in rejecting 5 per cent of the total number of potential emit- 
ters, but it should be noted that visual inspection of a sub-sample of these 
sources show that the majority are likely to be spurious, due to the combina- 
tion of a faint detection in NB921 and a < 3 a limit in z. Furthermore, the 
results in this paper remain unchanged even if these sources are included in 
the analysis. 



and 7 x 10~ 18 erg s -1 cm -2 for the NB921 data. The first-pass 
NBh sample of potential emitters contains 439 excess candidates 
out of all 23394 NB H detections in the entire NB H area (0.78 deg 2 ), 
while the NB921 sample has 8865 potential emitters out of 347341 
NB921 individual detections over the entire NB921 area coverage. 
Table |2]provides a summary of the number of sources and emitters 
throughout the selection proccess. 



3.2. 1 Visual inspection and star exclusion 

All NBh potential emitters are visually inspected in both the broad- 
band and narrow-band imaging. Twenty six (26) sources were re- 
moved from the sample as they were flagged as likely spurious. The 
majority of these (15) correspond to artifacts caused by bright stars 
that are on the edges of two or more frames simultaneously. The 
remaining 1 1 sources removed were low S/N detections in noisy 
regions of the NBh image. After this visual check, the sample of 
potential NBh is reduced to 413. 

Even by applying a conservative EW cut, the sample of poten- 
tial emitters can be contaminated by stars. Fortunately, stars can be 
easily identified by using the high-quality, multi-wavelength colour 
information available for the SXDF-UDS field. In particular, an op- 
tical colour vs. near-infrared colour (e.g. B — z vs. H — K) is 
able to easily separate stars from galaxies. Here, sources satisfy- 
ing (B - z) > 13.5(H - K) + 2.0 (AB) are classed as stars. 
By doing this, 2 potential stars are identified in the sample of po- 
tential emitters selected from the NBh data, and 118 (out of 5623 
emitters) from the NB921 data (the higher number of stars in the 
NB921 sample is driven both by a much larger sample size and, 
more importantly, because of the lower EW cut used to probe down 
to weak [Oil] EWs). These are excluded from the following analy- 
sis. It should be noted that all these sources present (B — z) colours 
significantly larger than 13.5(if — K) + 2.0 (AB), and therefore 
the identification of these as potential stars is not affected by small 
changes in the separation criteria. One of the potential NBh emit- 
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Figure 4. Left: photometric redshift distribution (peak of the probability distribution function for each source) for the NBh emitter candidates. Right: the 
equivalent distribution for NB921 emitter candidates. Both distributions peak at z ~ 1.47, corresponding to Ha: and [On] emitters, respectively, but other 
populations of emitters are also found, such as Ha at z — 0.4 and H/3/[Om] at z ~ 0.83 for the NB921 data and Paschen-lines at z < 1, [Olll] at z ~ 2.2 
and [On] at z ~ 3 for the NBh data. Lycn emitter candidates in the NB921 dataset (51 with reliable photometric redshift within the SpUDS coverage) with 
photo-z of z ~ 6 — 7 are not shown. 



ters flagged as a star is likely to be a cool T-dwarf (see Sobr al et al.| 
|2009b) . 

The final sample of NBh emitters over the entire NBh area 
contains 411 sources, of which 135, 69, 136 and 71 are found over 
the NE, NW, SE and SW fields, respectively (« 530 deg" 2 over the 
entire field down to the average NBh UDS depth). Note that by re- 
stricting the analysis to the UKIDSS-UDS area coverage (matched 
to the NBh and NB921 simultaneous coverage), which will be 
mostly used throughout this paper, the survey covers 0.67 deg 2 ; the 
sample of NBh emitters is reduced to 295 sources (the rest of the 
sources are outside this area), while the sample of NB921 emitters 
has 5505 sources. 



3.3 Distinguishing between different line emitters 

There are many possible emission-lines which can be detected in- 
dividually by the NBh and the NB921 filters. For z < 1 galaxies, 
the NBh filter is sensitive to lines such as Pa/3 at z = 0.26 or Pa7 
and Hel at z = 0.49 or Hell at z = 0.96, while at z > 1 the 
(main) possible emission lines are Ha at z ~ 1.47, [Olll]/H/3 at 
z ~ 2.25 and [Oil] stz = 3.34, among others. The NB921 filter is 
mostly sensitive to Ha at z ~ 0.4, [Olll]/H/3 at z ~ 0.83, [Oil] at 
z ~ 1.47 and hya at z ~ 6.6. Ha-[Oll] line emitters at z = 1.47 
can be selected with a significant narrow-band excess in both bands, 
and thus simultaneous excess sources provide an extremely robust 
means of selecting a z = 1.47 sample. This is explained in §3.4| 
First, in order to select [Oil] emitters below the Ha flux limit (see 
Figure [2|, to assess the robustness and contamination of different 
selections, to investigate the range of different emitters, and to al- 
low the selection in sky areas where coverage is not available in 
both filters, alternative approaches are considered. 




0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 
Spectroscopic Redshift 

Figure 5. A comparison between the best photometric redshifts in the UDS 
field and the robust spectroscopic redshift for such sources in UDS with 
spectroscopic redshift. NBh excess sources with spectroscopic redshifts 
are all identified with a known emission line falling into the NBh filter. 
It is noteworthy that the limited (and very biased towards the strongest line 
emitters and AGN) spectroscopic sample reveals that some confirmed Ho; 
emitters have photo-zs which place them at a higher and lower redshift - 
the use of NB921 data is able to identify these emitters which the photo-zs 
could miss. Vertical lines indicate the redshifts of the main emission lines 
detected by the NB H filter, specifically [Nl] at z = 0.545, [Svm]9914 
at z = 0.625, [On] 7625 at z = 1.125, Ha at z = 1.47, [Olll] 5007 at 
z = 2.23 and H/3 at z = 2.33. 



3.3.1 Photometric redshift analysis 

Multi-wavelength data can be used to effectively distinguish be- 
tween line-emitters at different redshifts (different emission lines) 
and separately obtain Ha and [Oil] selected samples of galaxies at 



z = 1.47 from the NB H and NB921 data-sets. Out of 5505 NB921 
potential emitters (excluding stars) within the UKIDSS UDS area, 
2715 (« 50%) have a > 3a detection in the near-IR UKIDSS 
data (K < 25 AB), while all NB H emitters are detected in K. 
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Figure 6. Top: The i— K vs. B— (AB) colour-colour distribution of 
NB921 emitters, clearly separating the 3 types of emitters that make the 
bulk of the sample of emitters. Simple lines for distinguishing these emit- 
ters - supported by evolution tracks (see S09) - are also shown. Bottom: 
An equivalent plot for the NBh emitters (AB), again demonstrating that 
lower redshift emitters can be easily isolated from the z = 1.47 Ho; emit- 
ters. Higher redshift emitters, however, occupy a similar region, and can not 
be robustly isolated with this set of colours alone. The two pivot points for 
the z ~ 0.4 and z ~ 0.8 separation are [0.63,0.84] and [1.71,2.77]; while 
the two pivot points used to separate z ~ 0.8 from z <~ 1.5 are [0.26,0.89] 
and [1.17,3.39]. 

Robust photometric redshift^] are only available over a matched 
0.64 deg 2 area, mostly due to the overlap with the Subaru data and 
the deep Spitzer coverage (SpUDS, PI J. Dunlop), and for sources 
with K < 23 (M. Cirasuolo et al., in prep.). Photo-zs are derived 
using a range of photometry (UBVRizJHK , together with IRAC 
3.6 and 4.5 bands) and are available for 25^] NBh potential line 
emitters and 1021 NB921 excess sources. 



5 The photometric redshifts for UDS present a(Az) = 0.015, where Az 
= ( z phot — Zspec)/ (1 + Zspec)- The fraction of outliers, defined as sources 
with Az > 3<j (Az), is lower than 2%. These photometric redshifts were 
used in |Sobral et al^(2010|[20rT) . 

6 The area reduction by itsell results in the loss of 297 NB921 emitters 
and reduces the sample of NBh line emitters from 41 1 to 295. It should be 
noted, nonetheless, that faint sources (with no reliable photo-z information) 
and those outside the photo-z area can still be further investigated by using 
colour-colour criteria. 



Figure|4]shows the photometric redshift distribution for the se- 
lected narrow-band emitters in NBh (left panel) and NB921 (right 
panel) with available photometric redshifts in UDS, demonstrating 
the common peak at z ~ 1.4 — 1.5, associated with the Ha/[Oll] 
lines being detected in each narrow-band filter. In addition to this, 
the other peaks are easily identified as Ha at z = 0.4 and [Olll]/H/3 
at z = 0.83 for the NB921 emitters. 



3.3.2 Spectroscopic redshift analysis 

Although an extensive spectroscopic sample is not yet available in 
the UKIDSS-UDS field (O. Almaini et al. in prep.; H. Pearce et 
al. in prep.), matching the samples of emitters with all the spec- 
troscopic redshifts published in the literature in this field ( Yamada 
et al.|2005(|Simpson et al.|2006[|Geach et al.|2008Uvan Breuke 



len et al. 2007, Ouchi et al. 2008 |Smail et al.||20081 |Ono et al 



2010^ allows the spectroscopic confirmation of 17 Ha emitters at 
z = 1.455 - 1.48 (over the full 0.78 deg 2 area), two [Olll]5007 
emitters at z = 2.23, an H/3 emitter at z = 2.33, and three lower 
redshift emitters ([Nl] at z = 0.545, [Svill] 9914 at z = 0.625 and 
[Oil] 7625 at z = 1.125), for the NB H data (see Figure[5]). Further- 
more, follow-up spectroscopy with SINFONI on the VLT of 6 Ha 
candidates has resulted in confirming all those sources, with spec- 
troscopic redshifts 1.45 — 1.47 (A. M. Swinbank et al. in prep.), 
and observations with FMOS on Subaru have confirmed a further 8 
Ha emitters (E. Curtis-Lake et al., in prep.), resulting in a total of 
31 spectroscopically confirmed z ~ 1.47 Ha emitters. 

For NB921 emitters, the limited spectroscopy confirms nine 
z = 0.4 Ha emitters, 12 z = 0.83 [Om] 5007 emitters, a z = 1.1 
H7 emitter, ten 4000A breaks (in which the NB921 filter probed 
light just to the red of the break, while the z band is dominated by 
emission blueward of the break, resulting in an excess in NB921), 
eight [Oil] emitters at z = 1.47, nine z = 2.31 Mgll emitters 
(AGN) and a z = 3.87 CIII] emitter, among others. 



3.3.3 Colour-colour separation of emitters 

Colour-colour diagnoses can be valuable tools to explore the ex- 
tremely deep broad-band photometry available, particularly for 
sources with no photometric redshift information (either because 
they are faint in the if -band, or because they are found outside the 
matched Subaru-Spitzer area). In S09, the BRiK colour-colour di- 
agram (B — R vs. % — K) is shown to isolate z ~ 0.8 emitters from 
lower and higher redshift emitters. As Figure [6] shows, the same 
colour-colour separation is also suited to separate z ~ 1.5 emitters 
from lower redshift emitters. Such colour-colour space is found to 
be particularly suited to distinguish between the bulk of the NB921 
emitters, as Figure [6] shows, clearly isolating Ha, [Olll]/H/3 and 
[Oil] emitters. Based on the spectroscopic confirmations and the re- 
liable photometric redshifts, empirical colour-colour selection cri- 
teria are defined (see Figure[6]) to distinguish between emitters. The 
same separation line produces relatively clean samples of z > 1.3 
emitters for the NBh sample. However, for this filter the Ha emit- 
ters at z rsj 1.5 and the [Olll]/H/3 emitters at z ~ 2 — 3 have similar 
BRiK colour-colour distributions, and a new set of colours needs 
to be explored to separate z ^ 1.5 and z > 2 emitters, after using 
the BRiK technique to remove low-z emitters. As Figure[7]shows, 
i — z vs. z — K colours provide a good separation between z ~ 1.5 
and z > 2 emitters (the selection is similar to the widely used BzK 



See UKIDSS UDS website for a redshift compilation by O. Almaini. 
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Figure 7. The z— K vs. i— z (AB) colour-colour distribution of NBh 
emitters, demonstrating a method to distinguish between the population of 
z ~ 1.47 Hen emitters, which dominate the sample, and higher redshift 
emitters such as [Om] and H/3 emitters at z ~ 2.2. The dashed line presents 
the simple separation used (in combination with the removal of lower red- 
shift emitters using the criteria in Figure[6) for sources with no photometric 
redshifts. The colour separation is similar to the BzK selection, but by using 
the i — z colour it better picks up the 4000A break for z = 1.47 galaxies, 
thus separating them better from z ~ 2 galaxies than using B — z. The two 
pivot points for the z ~ 2.2 and z ~ 1.47 separation are [0.09,1.50] and 
[0.49,3.96]. 



method, but can separate z ~ 1.5 galaxies from those at 2 ~ 2 and 
higher much more effectively), defining a colour-colour sub- space 
where Ha emitters at z = 1.47 should be found. 

These results demonstrate that a relatively good selection of 
z rsj 1.5 emitters can be obtained using just 5-band (BRizK) 
photometry. The BRiK selection can be used to remove low-z 
emitters, and a further izK analysis is capable of removing higher 
redshift emitters (although there's little contamination of the latter 
in the NB921 sample). 
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Figure 8. The photometric redshift distribution of simultaneous NBh & 
NB921 narrow-band emitters when compared to all NBh emitters within 
the NB921 imaging area. By matching both samples, it is possible to define 
a robust sample of Hai&[Oii] emitters without the need for colour or photo- 
metric redshift selection. Note that emitters with K < 23 AB are not found 
in the photo-z catalogue, and therefore not shown in the figure. 



Table 2. A summary of the number of sources within the different samples, 
including detections, candidate emitters, those within the matched area, 
those for which the full set of colours are available, sources with available 
photometric redshifts, and those selected as z = 1.47 emitters. 



Sample 


NB H 


NB921 


NB Detections (> 3 a, full area) 


23394 


347341 


Candidate emitters (full area) 


439 


8865 


After visual+star rejection (full area) 


411 


8747 


Matched area NB H +NB921 (0.67 sq deg) 


297 


5505 


Selection colours available (0.67 sq deg) 


297 


2715 


Photometric redshifts available 


257 


1021 


Selected z = 1.47 emitters (0.67 sq deg) 


190 (Ha) 


1379 ([Oil]) 



3.4 The dual narrow-band approach at z = 1.47 

As Figure [T] demonstrates, the two narrow-band filter profiles are 
extremely well matched in redshift when considering the detec- 
tion of the Ha and [Oil] emission lines (although the NB921 filter 
probes a slightly wider range in redshift). The match can be fully 
explored to select robust Ha emitters at z = 1.47, as these should 
be detected as [Oil] emitters in NB921 (because of the depth of 
those data). Within the matched NBh-NB921 area (out of a sample 
of 297 NBh emitters), a sample of 178 dual-emitter sources is re- 
covered. Figure[8]shows the photometric redshift distribution of the 
double emitters, compared with the distribution of all NBh emitters 
with photo-zs (for roughly the same area). The results clearly sug- 
gest that the dual-emitter candidate criteria is able to recover essen- 
tially the entire population of z = 1.47 Ha emitters. Moreover, the 
dual-emitter selection also recovers some sources with photometric 
redshifts which are significantly higher (3 sources) and lower (14 
sources) than z — 1.47 (accounting for the 1 a error in the photo- 
zs) and that would have been missed by a simple photo-z selec- 
tion. There are spectroscopic redshifts available for 7 out of the 14 
lower-photo-z sources, and they all confirm those sources to be at 
z = 1.45 — 1.48, indicating that their selection by the dual-emitter 
technique is very reliable. 



The sources for which the photo-z fails have high Ha fluxes, 
and/or are very blue. It is therefore likely that the photo-zs fail 
for these sources because they do not present identifiable breaks 
(pushing photo-zs to a lower redshift solution), and/or because they 
present strong enough Ha lines to contaminate the photometry. 



3.4. 1 Searching for low [ Oil ] EWs Ha emitters 

Figure |9] shows the distribution of [Oll]/Ha line ratios for the emit- 
ters selected in both narrow-band filters. In order to ensure maxi- 
mum completeness at the faintest [Oil] fluxes, and minimise any 
possible biases, a search for extra Ha emitters with significant 
(S > 3) NB921 colour-excess but low [Oil] EWs (i.e. waiving 
the 25 A EW requirement), is conducted. This yields 12 additional 
sources above the NB921 flux limit. All of these additional sources 
present colours and photometric redshifts consistent with being 
genuine z ~ 1.47 sources and probe down to the lowest [Oll]/Ha 
line fractions (see Figure[9]); their colours are consistent with being 
sources affected by higher dust extinction. These 12 sources are 
added to the robust double-emitter sample. 
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Figure 9. The distribution of the [On]/Hai line ratios for the robust sample 
of 190 simultaneous Hai-[Oii] line emitters at z = 1.47 (the shaded region 
represents the 12 sources with the low [On] EWs). The vertical dashed line 
represents the lowest line ratio that the combined survey is able to probe 
at the lowest Hen fluxes; this shows that the [On] survey is deep enough to 
recover essentially all Ho; emitters and that the line flux distribution is not 
a result of a bias caused by insufficient sensitivity at the lowest line flux 
ratios. 



3.4.2 The completeness of the Ha-[ Oil ] double selection 

The completeness limit on the [Oll]/Ha flux ratio (see Figure[9} in- 
dicates that the deep NB921 flux limit guarantees a very high com- 
pleteness of the dual-emitter sample, and suggests that the small 
number of sources that photo-zs indicate as being at z ~ 1.47, but 
which are not NB921 excess sources, are actually at other redshifts. 
To test this further, the photometric redshift and colour-colour se- 
lections are applied in the same matched sky area. These identify 
32 Ha candidate sources which are not selected by the dual-emitter 
approach. Four of these sources have spectroscopic redshifts avail- 
able, and all four of these are indeed spectroscopically confirmed 
to lie at different redshifts (cf. Figure [5]). Of the remaining 28, 18 
sources have significantly negative excesses in NB921, and the re- 
maining 10 have < 3S excesses, which would correspond to a 
range of NB921 line fluxes from - 0.1 - 1 x 10" 18 erg s" 1 cm" 2 . 
If real, these line fluxes would imply very low [Oll]/Ha flux ratios, 
indicating that the sources must be highly extinguished (Ah« ~ 
2 — 4 mag); however, their UV and optical colours appear incon- 
sistent with such high extinction (see §5.5). Thus, these sources are 
likely to be a mixture of other line emitters of weaker emission lines 
close to Ha (e.g. [Nil], [Sll]), which the photometric redshifts are 
not sufficiently accurate to distinguish (see Sobral et al. 2009a). 



3.5 Selecting robust Ha and [Oil] z = 1.47 emitters 

3.5.1 Ha emitters at z — 1.47 

§|34]shows that with the flux limits of this study, the NB H -NB921 
is clearly a clean and highly complete means of selecting z — 1.47 
Ha emitters. The sample of robust Ha emitters at z = 1.47 is 
therefore derived solely by using the dual-emitter selection, within 
the matched area. This results in the selection of very robust 190 
Ha- [Oil] sources at z — 1.47. 



3.5.2 [Oil] emitters at z = 1.47 

In order to select [Oil] emitters, in addition to the dual-emitter se- 
lection at bright fluxes, both photo-zs and(or) colours are used. In 
particular, the z = 1.47 [Oil] candidates are defined to be narrow- 
band excess sources with z m i n < 1.47 < z max (where z m i n and 
Zmax are the 3-a redshift limits of the principle peak in the photo- 
metric redshift probability distribution - photo-zs of emitters have 
a typical Az of « 0.13) or those that are found within the defined 
BR-iK colour-colour space ((B-R) < 0.26 or (i-K) > 4.92 
or (i - K) > 4.121(B - R) + 1.349; see Figure^. Even though 
the contamination from higher redshift (z > 1.5) emitters (such 
as Mgll) is likely to be low, the colour-colour selection presented 
in Figure|7]is also applied ((i — z) < 0.24 or (z — K) > 5.33 or 
(z — K) > 6.15(i — z) + l. 394), to remove potential higher redshift 
(z > 1.5) emitters. This results in a sample of 1379 [Oil] selected 
emitters down to the flux limit of the survey. However, as noted be- 
fore, many faint NB921 emitters are not detected in at least one of 
the bands necessary for the colour-colour identification, and there- 
fore are not included directly in the [Oil] sample. j ]4. 3 [ investigates 
how this can bias the determination of the faint-end of the lumi- 
nosity function and derives a correction to account for the sources 
which are missed. 



4 Ha AND [Oil] LUMINOSITY FUNCTIONS 

4.1 Contribution from adjacent lines 

While the NB921 filter can measure the [Oil] emission line (which, 
in fact, is a doublet) without contamination from any other nearby 
emission line, this is not the case for NBh and the Ha line. The 
nearby [Nil] emission lines can contribute to the NBh emission 
flux, and therefore both EWs and fluxes will be a sum of [Nil] 
and Ha. The limited spectroscopic follow-up in the H band with 
SINFONI (A. M. Swinbank et al. in prep.) provides good enough 
S/N and spectral resolution to compute individual [Nil] frac- 
tions/corrections. The data show a range of [Nil] fractions between 
0.1-0.4 and, despite being a limited sample, the results are broadly 
consistent with Equation 3 of S09, which estimates the [Nil] con- 
tamination as a function of total rest-frame EWo([Nll]+Ha), based 
on a large SDSS sample. For the analysis presented in this paper, 
the polynomial approximation used in S09 has been re-computed 
using a higher order polynomial which is able to reproduce the 
full SDSS relation between the average, log([Nll]/Ha), /, and 
log[EW ([Nll]+Ha)], E: f = -0.924 + 4.802£ - 8.892£ 2 + 
6.701£ 3 - 2.27E 4 + 0.279£ 5 . This relation is used to correct all 
Ha fluxe^Jin this paper. The average correction for the z = 1.47 
sample is 0.22 (average rest-frame EW (Ha+[Nll]) of 130A). 

4.2 Individual line completeness 

It is fundamental to understand how complete the samples are as a 
function of line flux. This is done using simulations, as described in 
S09. Briefly, the simulations consider two major components driv- 
ing the incompleteness fraction: the detection completeness (which 
depends on the actual imaging depth and the apertures used) and the 
incompleteness resulting from the selection (both EW and colour 
significance). The first component is studied individually per frame, 

8 The relation is used to re-compute Hen fluxes in S09 for full consistency; 
the average correction at z = 0.84 is 0.28 (which compares to 0.25 in S09). 
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Figure 10. A study of the completeness fraction (defined as the fraction of 
sources with a given flux recovered by the selection against the actual num- 
ber of sources with that flux) in order to address incompleteness as a func- 
tion of line flux for both NB H and NB921 samples in the SXDF-UKIDSS- 
UDS field. Note that both image detection and colour-magnitude selection 
are taken into account as these are both sources of incompleteness. This 
confirms the flux limit computed for both bands (down to ~ 3 a). Note 
that the NBh curve shown is used to correct the Hen luminosity function di- 
rectly for incompleteness (for each field), while the NB921-[On] function 
plotted is the one that is used when correcting the [On] luminosity func- 
tion (obtained by including both the flux incompleteness and the selection 
incompleteness - see Figure 1 1 and Section 4.3). 
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Figure 11. Fraction of [On] and non-[On] NB921 emitters (amongst those 
with K detections, allowing for a colour-colour separation) as a function of 
emission line flux. A significant number of potential line emitters are not in 
the final [On] colour-selected sample as they are not detected in the K band 
imaging; the fraction of all emitters with a K detection is also shown. 



20% of the size of the correction is added in quadrature to the 
other uncertainties to account for possible inaccuracies in the sim- 
ulations. 



by adding a set of fake galaxies with a given input magnitude to 
each frame and obtaining both the recovery fraction and the recov- 
ered magnitude. The second component is studied by using sources 
which have not been selected as emitters, and adding an emission 
line with a given flux to all those in order to study the fraction re- 
covered. In an improvement to S09, sources classed as stars and 
those occupying the z < 1 region of the BRiK diagram are not 
used in this simulation - this results in a set of galaxies with an 
input magnitude distribution which is very well-matched to the 
z — 1.47 population, providing a more realistic sample to study the 
line completeness of the survey for z — 1.47 Ha emitters (see Ap- 
pendix A, which quantifies the difference of using this methodfj 

For each recovered source, a detection completeness is asso- 
ciated, based on its new magnitude. The results (survey average) 
can be found in Figure [To] but it should be noted that because of 
the differences in depth, simulations are conducted for each indi- 
vidual frame, and the appropriate completeness corrections applied 
accordingly when computing the luminosity function. 

A similar procedure is followed for the NB921 data (average 
results are also shown in Figure [ToJ; corrections are also derived 
and used individually for each field, although the differences in 
depth are not as significant as for the NBh data. Similarly to the 
NBh analysis, non-emitter sources which are classed as stars and 
low redshift (z < 1) are rejected when studying how the complete- 
ness of the survey varies with line flux in order to better estimate 
how complete the survey is specifically for z = 1.47 line emitters. 

For any completeness correction applied, an uncertainty of 



9 The simulations done in S09 are repeated following this approach, and 
the results are used to re-compute the Hen luminosity function at z = 0.84 
in order to guarantee a completely consistent comparison. 



4.3 [On] selection completeness 

A significant fraction of NB921 -selected emitters do not have pho- 
tometric redshifts (79%), and half of the sample is undetected in K 
above the 3a limit. Therefore, these cannot be robustly classified, 
and they are not included in the [Oil] sample, even though at least 
some of them may well be genuine [Oil] emitters at z = 1.47. This 
can cause a potentially significant incompleteness in the [Oil] sam- 
ple, resulting in an underestimation of the luminosity function. In 
order to investigate this source of incompleteness and its effect in 
the determination of the luminosity function, the fraction of emit- 
ters with K detections (which generally sets the limit for colour- 
colour selection, as the other bands used are significantly deeper 
- though even in the optical bands typically 33% of the sources 
are undetected) is studied as a function of line flux. As Figure [TT] 
shows, the fraction of NB921 excess sources with K detections 
falls as a function of decreasing line fluxes, suggesting that the in- 
completeness is higher at the lowest fluxes. This can have impli- 
cations for determinations of the faint end slope of the luminosity 
function. 

Figure[TT]also shows that, for the classifiable sources, the frac- 
tion of [Oil] emitters rises with decreasing flux at the highest fluxes, 
and then seems to remain relatively flat at ~ 60% down to the 
lowest fluxes. Assuming that this same distribution is true for the 
unclassified sources (which is likely, since the [Oil] fraction also 
remains roughly constant with if -band magnitude for the faintest 
K magnitudes probed), it is possible to estimate a correction at a 
given flux that is given by the ratio of the classifiable sources to 
all sources at that flux. This is equivalent to including unclassified 
emitters in the luminosity function calculation with a weight that is 
given by the fraction of [Oil] emitters at that flux (see Figure 10). 
An uncertainty corresponding to 20% of this correction is added in 
quadrature to the other uncertainties. 
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4.4 Volume: Ha and [Oil] surveys 

Assuming the top-hat (TH) model for the NB H filter (FWHM of 
211.1 A, with A™ 1.606 um and A^ax 1.627 nm), the 
Ha survey probes a (co-moving) volume of 2.667 x 10 5 Mpc 3 . 
When matched to the NB921 coverage, this is reduced to 2.2872 x 
10 5 Mpc 3 due to the reduction in area - the matched FWHM is the 
same. While this is the total volume probed for the common flux 
limit over the entire coverage, each pointing reaches a slightly dif- 
ferent flux limit, and therefore at the faintest fluxes the volume is 
smaller (as only one pointing is able to probe those). This is taken 
into account when determining the luminosity function: only areas 
with a flux limit above the flux (luminosity) bin being calculated 
are actually taken into account. In practice, this results in only us- 
ing no one quarter of the total area for the faintest bin (with this 
being derived from the deeper WFCAM pointing). 

Assuming a top-hat model for the NB921 filter, (FHWM 
of 132 A, with A™ 0.9130 [im and A™ x 0.9262 ^im), 
the [Oil] survey probes a volume (over 0.67 deg 2 ) of 2.5102 x 
10 5 Mpc 3 when assuming a single line at 3 727 A, and a volume 
of 2.6363 x 10 5 Mpc 3 using the fact that the [Oil] line is actually a 
doublet; 3726.1 A and 3728. 8A. For simplicity, because the change 
in volume is less than 5 per cent (and therefore much less than the 
errors), and for consistency with other authors (allowing a better 
comparison), the volume used assumes a single line. 



4.5 Filter Profiles: volume and line ratio biases 

Neither NBh or NB921 filter profiles are perfect top-hats (see Fig- 
ure [TJ. In order to evaluate the effect of this bias on estimating the 
volume (luminous emitters will be detectable over larger volumes - 
although, if seen in the filter wings, they will be detected as fainter 
emitters), a series of simulations is done. Briefly, a top-hat volume 
selection is used to compute a first-pass (input) luminosity function 
and derive the best Schechter Function fit. The fit is used to gener- 
ate a population of simulated Ha emitters (assuming they are dis- 
tributed uniformly across 1.40 < z < 1.52); these are then folded 
through the true filter profile, from which a recovered luminosity 
function is determined. Studying the difference between the input 
and recovered luminosity functions shows that the number of bright 
emitters is underestimated, while faint emitters are slightly overes- 
timated (c.f. S09 for details). This allows correction factors to be 
estimated - these are then used to obtain the corrected luminosity 
function. 

Figure [T] also shows how the NBh and NB921 filters are very 
well matched, although the [Oil] coverage is slightly wider than 
the Ha coverage. In order to evaluate how this might affect the 
results on line ratios, a series of simulations is done. Simulated 
[Oll]+Ha emitters are distributed uniformly in a volume defined 
by 1.40 < z < 1.52, which contains the entire transmission re- 
gions of both profiles. Emitters are given a wide range of Ha fluxes 
based on the observed Ha luminosity function, and [Oil] fluxes 
corresponding to [Oll]/Ha line ratios between and 2.0. The real 
filter profiles are then used to recover, for each emitter, both the 
[Oil] and Ha fluxes, and therefore allow the study of both the re- 
covered line fluxes and the recovered line ratios. Based on these 
results - presented in Figure [12]- the line ratios should be accurate 
within ~ 20%. 
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Figure 12. Top: the variation of the measured [On]/Hai line ratio as a func- 
tion of galaxy redshift due to the small differences in the transmission func- 
tion of both filter profiles. Note that the flux limit (in units of erg s — 1 cm -2 ) 
also varies very significantly as a function of redshift, and thus line ratios 
are likely to be correct to within less than 20% for the large bulk of the sam- 
ple. Bottom: by using the Ho; luminosity function to produce a population 
of galaxies distributed over a wider redshift range than the filter profile, it is 
possible to use both NBh and NB921 filter profiles to study in a simple way 
the distribution of the ratio between the recovered and input line fractions. 
The results recover the input line fraction with 17% standard deviation and 
with less than 1% recovered with a line fraction increased by more than 
+50%. 



4.6 Extinction Correction 

The Ha emission line is not immune to dust extinction, although 
it is considerably less affected than the [Oil] emission line. Mea- 
suring the extinction for each source can in principle be done by 
several methods, ranging from spectroscopic analysis of B aimer 
decrements to a comparison between Ha and far-infrared deter- 
mined SFRs, but such data are currently not available. 

In this Section, the analysis of the Ha luminosity function is 
done assuming Ah« = 1 mag of extinction at Hap°) as this allows 
an easy comparison with the bulk of other studies which have used 
the same approach. A few studies have used a Ha dependent ex- 
tinction correction - either derived from Garn et al. ( 2010 ), or from 
Hop kins et al.| ( |200T] ). However, ^suggests that at least the over- 
all normalization of such relation at z ~ 1.5 is significantly lower 
(dust extinction is not as high as predicted for the very high lumi- 
nosities probed) and therefore over-predicts dust-extinction correc- 
tions by ~ 0.5 — 0.7 mag (overcorrecting luminosities by 0.2 — 0.3 
in dex). Finally, the [Oil] luminosity function is presented without 
any correction for extinction, in order to directly compare it with 
the bulk of other studies. Detailed extinction corrections and a dis- 
cussion regarding those are presented in Section[5] 



10 Corresponding to ~ 1.76 mag of extinction, or a factor of 5x at [On] 
for a Calzetti et al. (2000) extintion law. 
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Figure 13. The Hen luminosity function evolution up to z ~ 2.3. The new 
measurement at z = 1.47 is combined with HiZELS measurements (S09 
and G08) and other narrow-band stud ies (Gallego et al.|1995[|Ly et al.|2007[ 
IShioya et al.|2008l|Geach et al.|2008||Hayes et al.|2010||Tadaki et aLpoTTj . 
This confirms the strong evolution from z = 0toz = 2.23. 
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Figure 14. The [On] luminosity function derived from this study, at z = 
1.47 and a comparison with measurements at different redshifts from Gal- 
|lego et al.| ( [2002) ; |Takahashi et al.| ( |2007) ; |Ly et al.|j2007> ; |BayEss et al.| 
(2011) and with a measurement at the same redshift (using the same narrow- 
band filter, but on a different field) by Ly et al. (2007). There is a clear 
evolution in the [On] luminosity function, at least up to z ~ 2. 



4.7 Ha Luminosity Function at z = 1.47 and Evolution 

By taking all Ha selected emitters, the luminosity function is com- 
puted. As previously described, it is firstly assumed that the NBh 
filter is a perfect top-hat, but the method of S09 is applied to cor- 
rect for the real profile (see Section [43| ). Candidate Ha emitters are 
assumed to be all at z = 1.47 (as far as luminosity distance is con- 
cerned). Results can be found in Figure [13] Errors are Poissonian, 
with a further 20% of completeness corrections added in quadra- 
ture. The luminosity function is fitted with a Schechter function 
defined by three parameters a, 0* and L*: 

0(L)dL = <A*(^)V^*>d(A). (5 ) 

In the log forrrplthe Schechter function is given by: 

</>(L)dL = In 10 0* (-^y e- (L/L ' ] (-^) dlogL. (6) 

A Schechter function is fitted to the Ha luminosity function, with 
the best fit resulting in: 

log L^ a = 42.50 ± 0.23 erg s _1 
log<& a = -2.44 ± 0.33 Mpc -3 
a Hc , = -1.6 ±0.4 

The best-fit function is also shown i n Figure |T3] together with the 
luminosity functions determined by |Ly et al.||2007| > - extending 
the work by Gallego et al. (1995) - at z ~ and that from Shioya 
|et al.| ( [2008] ) at z — 0.24. Also presented on Figure[l3]are luminos- 
ity functions from the other HiZELS redshifts. The new luminos- 
ity function uses the revised catalogues presented in Sobral et al. 
(2010) and Sobral et al. (201 1), and uses new completeness correc- 
tions, re-computed fluxes and the rest-frame EW dependent [Nil] 
correction. The changes to the luminosity function are only minor. 
At z = 2.23, the luminosity function has also been re-calculated 
by combining the results from Geach et al. (2008), with new ultra- 
deep measurements of the faint end by |Hayes et al.| ( |2010| ) - using 

1 1 log 10 ; note that the definition of </>* used in S09 and Geach et al. (2008), 
4>* s , is given by </>*=</>* x In 10. 



HAWK-I on VLT - and Tada ki et al.|{20ll) - using MOIRCS on 
Subaru. Best fits for the HiZELS luminosity functions at z — 0.84, 
z — 1.47 and z — 2.23 are presented in Table|3] 

The results confirm the strong evolution from the local Uni- 
verse to z — 2.23 and provide further insight. In particular, while 
there is significant evolution up to z ~ 0.8, the Ha luminos- 
ity functions at z ~ 0.8, z ~ 1.47 and z ~ 2.23 (including 
the ultra-deep measurement by |Hayes et"aL]|2010| ) seem to agree 
well at the lowest luminosities probed (below logLna ~ 42.0) 
in normalisation and slope - all consistent with a relatively steep 
value of a ~ —1.6) - a similar a is found by |Ly et aT]p011| ). 
In contrast, at the bright end, L^ is clearly seen to continue to 
increase from z — 0.84 to z — 1.47 and to z = 2.23, with 
log L* oc 0.6z. This implies that the bulk of the evolution from 
z — 2.23 to z ~ 1 is happening for the most luminous Ha emitters 
(Lhq > 10 42 ergs -1 ), which greatly decrease their number den- 
sity as the Universe ages: the faint-end number densities seem to 
remain relatively unchanged. 



4.8 The [On] Luminosity function at z = 1.47 and Evolution 

The [Oil] luminosity function at z — 1.47 (not corrected for dust 
extinction) is shown in Figure^] Note that this includes correc- 
tions for i) incompleteness in both detection and flux selection, ii) 
incompleteness in the redshift selection of [Oil] emitters (see Fig- 
ure 10 for the combined correction) and iii) correction for the filter 
profile not being a perfect tophat. Errors are Poissonian and include 
(in quadrature) 20% of the completeness correction applied). The 
fully-corrected [Oil] luminosity function is found to be well-fitted 
by a Schechter function, with the best fit resulting in: 

logL* OII] =41.71 it 0.09 ergs -1 
log0* OII] = -2.01 ±0.10 Mpc -3 
oqoii] = -0.9 ±0.2 

The results are compared with other studies at different red- 
shifts ( |Hogg et al.|1998l|Gallego et al.|2002||T5tahashi et al.|2007| 
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Table 3. The luminosity function and star-formation rate density evolution at the peak of the star formation history as seen by HiZELS at z = 0.84, z = 1.47 
and z = 2.23; assuming 1 mag extinction at Ha. Columns present the redshift, break of the luminosity function, L* H(x , normalisation, 4>* Ha and faint-end 
slope of the luminosity function, a. The three right columns present the star formation rate density at each redshift based on integrating the luminosity function 
down to the given luminosity limit (in log erg s _1 ). Star formation rate densities include a correction for AGN contamination of 10% at z = 0.84 (c.f. Garn 
et al. 2010) and 15% at both z = 1.47 and z = 2.23. 
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z = 0.84 ±0.02 


42.25 ±0.16 


-2.36 ±0.29 


-1.66 ±0.34 


0.03 ±0.01 


0.12 ±0.04 


0.13 ±0.06 


z = 1.47 ±0.02 


42.50 ±0.23 


-2.44 ±0.33 


-1.6 ±0.4 


0.06 ±0.02 


0.15 ±0.06 


0.16 ±0.05 


z = 2.23 ±0.02 


43.07 ±0.24 


-2.93 ±0.41 


-1.60 ±0.21 


0.12 ±0.06 


0.20 ±0.09 


0.22 ±0.10 



|Ly et al.|2007l|Bayliss et al.|2 011) - all uncorrected for dust ex- 
tinction - and shown in Figure [14] The comparison reveals a sig- 
nificant evolution in the [Oil] luminosity function from z ~ to 
z — 1.47. Such evolution seems to be simply described by a 0* 
and L* evolution up to z ~ 0.9 and a continuing L* evolution from 
z rsj 0.9 up to z ~ 2, in line with the results for the evolution of the 
Ha luminosity function. It should also be noted that there is good 
agreement with the |Ly et a l. (2007) luminosity function at the same 
redshift derived from a smaller area, but at a similar depth. 

The correction for redshift/emitter selection is found to be par- 
ticularly important at faint fluxes, setting the slope of the faint-end 
of the luminosity function (although it has little effect on the val- 
ues of L* and 0*). If no correction for unclassified sources is ap- 
plied, the best-fit Schechter function yields a « —0.2 (compared 
to the corrected best- fit of a « —0.9). In the extreme case where all 
the non- selectable emitters are assumed to be [Oil], it would yield 
a « —1.8. This makes it clear that the crucial data that one needs 
in order to improve the determination of the faint-end slope is not 
new, deeper NB921 data, but rather significantly deeper B, i, z and 
K data, which will enable to completely distinguish between dif- 
ferent line emitters at faint fluxes (or at least provide a more robust 
correction by classifying a much higher fraction of emitters). 

4.9 Evolution of the Ha and [On] Luminosity Functions 

Robust measurements of the evolution of both the Ha and [Oil] 
luminosity functions up to z ~ 2 have been presented. The results 
reveal that there is a strong, and consistent, evolution of both lumi- 
nosity functions. Moreover, while the evolution up to z ~ 1 can 
be described as both an increase of the typical luminosity (L*) and 
an increase of the overall normalization of the luminosity functions 
(0*), at z > 1 the bulk of the evolution can be simply described 
as a continuous increase in L* given by log L* oc 0.6 z. The cur- 
rent results also show that the faint-end slope (and the number den- 
sity of the faintest star-forming galaxies probed) seems to remain 
relatively unchanged during the peak of the star formation history 
(z rsj 1 — 2). Of course, the latter does not imply, at all, that the 
faint population is not evolving. Furthermore, particularly beyond 
z ~ 1, AGN contamination at the highest luminosities could still be 
polluting the view of the evolution of the star-forming population. 

4.10 The star formation rate density at z = 1.47 

The best-fit Schechter function fit to the Ha luminosity function 
can be used to estimate the star formation rate density at z = 1.47. 
The standard calibration of Kennicutt (1998) is used to convert the 
extinction-corrected Ha luminosity to a star formation rate: 

SFR(M yr -1 ) = 7.9 x 10~ 42 L Hc * (ergs -1 ), (7) 



which assumes continuous star formation, Case B recombination 
at T e = 10 4 K and a Salpeter initial mass function ranging from 
0. 1-100 Mq. A constant 1 magnitude of extinction at Ha is as- 
sumed for the analysis, as this is a commonly adopted approach 
which is a good approximation at least for the local Universe. ^5] 
will investigate this further based on the analysis of emission-line 
ratios, revealing that Ah« ~ 1 mag for the sample presented in 
this paper. A 15% AGN contamination is also assumed (Garn et 
al. 2010 found AGN ~ 10% contamination at z = 0.84, but 
the contamination is likely to be higher at higher redshift and 
the slightly higher flux limit of this sample). Down to the survey 
limit (after extinction correction, Lh« = 10 42 ergs -1 ), one finds 
Psfr = 0.07 =b 0.01 M yr -1 Mpc -3 ; an integration extrapolat- 
ing to zero luminosity yields /?sfr = 0.16=b0.05 M yr -1 Mpc -3 . 
Furthermore, by using the standard Kennicutt (1998) calibration of 
[Oil] as a star formation tracer (calibrated using Ha): 

SFR(M year -1 ) = 1.4 x 10" 41 L [0 n] (ergss -1 ), (8) 

it is also possible to derive an [Oil] estimate of the star forma- 
tion rate density at the same redshift, by using the complete in- 
tegral of the [Oil] luminosity function determined at z = 1.47, and 
assuming the same 1 magnitude of extinction at Ha. This yields 
Psfr = 0.17 =b O.O4M yr -1 Mpc -3 , in very good agreement 
with the measurement obtained from Ha. 

By taking advantage of the other HiZELS measurements, and 
other Ha based measurements, it is possible to construct a full and 
consistent view of the Ha-based star formation history of the Uni- 
verse. This is done by integrating all derived luminosity functions 
down to a common luminosity limit of Lhc*=10 41 5 erg s -1 . The re- 
sults are presented in Figure [15] revealing a strong rise in the star 
formation activity of the Universe up to z ~ 1 and a flattening or a 
small continuous increase beyond that out to z > 2. 



5 THE MATCHED Ha-[OII] VIEW 

This section presents a detailed comparison of [Oil] and Ha lu- 
minosities and line flux ratios (observed and corrected for dust- 
extinction) as a function of galaxy colour, mass and luminosity for 
the sample of matched Ha and [Oil] emitters at z = 1.47. These 
are compared with equivalent results for a similarly- selected sam- 
ple at z ~ 0.1 drawn from the SDSS. 

5.1 A SDSS sample at z - 0.1 

In order to compare the z = 1.47 Ha+[Oll] sample with a large 
local sample and provide a further insight into any important corre- 
lations between line fractions and dust-extinction, mass and colour, 
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Figure 15. The star formation rate density and its evolution with redshift up to z ~ 2.3 using Hen only, but also including the [On] measurement presented in 
this paper (assuming 1 magnitude of extinction at Hen) and compared to estimates at different redshifts from the literature (e.g. Hopkins 2004 Ly et al. 2011 
and references herein). Darker circles represent the results from Ho; studies which are also shown in Figure 13. This confirms a strong evolution in the star 
formation rate density over the last ~ 10 Gyrs, with a flattening or slow decline over 1 < z < 2 and a sharp decrease from z ~ 1 to z ~ 0. 



an SDSS sample was used. Data were extracted from the MPA 
SDSS derived data products catalogue^] The sample was defined 
to emulate a narrow-band slice at 0.07 < z < 0.1 (chosen to be 
distant enough so that the fibers capture the majority of the light 
- typically ^3 — 4 kpc diameter, but not too far out in redshift 
to guarantee that the sensitivity is still very high and the measure- 
ments are accurate). The sample was selected by further imposing 
a requirement that > 10 40-6 ergs -1 (observed, not aperture 
or dust extinction corrected); this guarantees high S/N line ratios 
and that the vast majority of sources are detected in [Oil] as well, 
allowing unbiased estimates of the line ratio distribution, as well 
as the detection of other emission lines (Olll], H/3 and [Nil]) that 
can be used to distinguish between AGN and star-forming galax- 
ies (c.f. |Rola et al.| 19971 |Brinchmann et al.|2004| ). Among 17354 
SDSS z ~ 0.1 Ha emitters, 498 were classified as AGN, implying 
a rsj 3 per cent contamination. Potential AGN were removed from 
the sample. Emission line fluxes are aperture corrected following 
a similar procedure as in Garn & Best (2010) - i.e., by using the 
ratio between the fiber estimated mass and the total mass of each 
galaxy. For 32 galaxies where the catalogued fiber mass is higher 
than the catalogued total mass, an aperture correction factor of 1.0 
was assigned. For those in which the fiber mass had not been deter- 
mined, the average correction for the total mass of that galaxy was 
assigned. The median fraction of flux within the SDSS apertures is 
32%. Note that the even though these aperture corrections change 
the total fluxes, line ratios remain unchanged. Finally, in order to 
provide a more direct comparison with the sample at z = 1.47, an 



EW cut of 20 A in Ha was applied (to mimic the selection done 
at z = 1.47), and galaxies with lower EWs were excluded (1701 
galaxies). The final SDSS sample contains 14451 star-forming Ha 
selected galaxies at 0.07 < z < 0.1. 

5.2 Calibrating [On]/Ha line ratio as a dust extinction probe 

Due to the difference in rest-frame wavelength of the [Oil] and Ha 
emission lines, and both emission-lines being tracers of recent star- 
formation, the [Oll]/Ha line ratio is sensitive to dust-extinction, 
even though metallicity can affect the line ratio as well. After cor- 
recting for dust extinction, several studies find a [Oll]/Ha average 
of « 1.0 - 1.4 in the local Universe (c.f. |Kewley et al.|2004) . 

The Ha/H/3 line ratio (B aimer decrement) is widely used as an 
extinction estimator, particularly up to z ~ 0.4, as it is relatively 
easy to obtain both emission lines. As the SDSS-derived sample 
is able to obtain reliable fluxes for the [Oil], H/3 and Ha emis- 
sion lines, it is possible to investigate (and potentially calibrate) 
[Oll]/Ha as a dust-extinction indicator, using the B aimer decre- 
ment directly. As Figure [T6] shows, [Oll]/Ha is relatively well cor- 
related with Ha/H/3, indicating that it is possible to use [Oll]/Ha to 
probe dust extinction within the observed scatter. For each galaxy 
in the SDSS-derived sample at z ~ 0.1, Ha/H/3 line fluxes are 
measured and used to estimate the extinction at Ha, Ah q , by us- 
ing: 



2.5/chc , / 2.86 \ 

: login I — t=tzz 1 



V Ha/H/3 J 



(9) 



See http://www.mpa-garching.mpg.ed/SDSS/ 



where 2.86 is the assumed intrinsic Ha/H/3 line flux ratio, appro- 
priate for Case B recombination, temperature of T = 10 4 K and 
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Figure 16. The variation of [On]/Ha line ratios as a function of Hck/H/3 
ratios for SDSS (balmer decrement), showing that they correlate well, and 
therefore it is possible to calibrate the observed [Oll]/Ha line ratio as a 
dust extinction indicator. The best polynomial fit to the data is also shown, 
together with the prediction from the Calzetti et al. extinction law. 



an electron density of n e — 10 2 cm -3 |Brocklehurst| 1 97 1 \ . The 
|Calzetti et al.| < |2000| > dust attenuation law is used to calculate the 
values of k\ = A\/ E{B — V) at the wavelengths of the Ha and 
H/3 emission lines, resulting in: 



A Hq = 6.531 log 10 Ha/H/3 - 2.981. 



(10) 



By using the Calzetti law, it is also possible to write a similar rela- 
tion by using [Oll]/Ha: 



A Hq = -2.2961og 10 [OII]/Ha + 2.291og 10 (R ), 



(11) 



where Ro is the unknown intrinsic [Oll]/Ha, but the slope of the 
relation is fully determined. Figure [16] shows that the Calzetti law 
matches the global trend well, but fails to predict the fine details. 
The failure to match the data perfectly results from a combination 
of different factors. The Calzetti et al. law does not include im- 
portant variations of the intrinsic [Oll]/Ha line fraction with e.g. 
metallicity. The Calzetti et al. extinction law is also based on a rel- 
atively small number of local galaxies, and, most of all, it is based 
on continuum light, and not on emission-lines. However, it is possi- 
ble to obtain a much better fit to the relation between [Oll]/Ha (\E0 
and Balmer decrement with a 4th order polynomial (shown in Fig- 
ure 16), and to derive an empirical relation between Ana directly 
from HafH/3: 

Arc = -4.30^ 4 - 11.30^ 3 - 7.39^ 2 - 2.94^ + 0.31 (12) 

with a scatter of 0.28 mag (in Ah q ). The scatter may be driven by 
the range of metallicities within the sample, as galaxies will have 
slightly different intrinsic [Oll]/Ha line ratios. 

The SDSS data can be used to investigate how the offset of a 
galaxy from the [Oll]/Ha vs Ha/H/3 relation depends upon metal- 
licity, and thus it is possible to measure the extent to which the 
[Oll]/Ha line ratio can be used to probe extinction, without metal- 
licity biases. Here, the 03N2 indicator ( jPettini & Pagel||2004| > is 
used as a tracer of metallicity (the gas-phase abundance of oxygen 
relative to hydrogen), computed by using: 

12 + log 10(O/H) = 8.73 - 0.32 x log 10 (O3H/3/N2Ha), (13) 
where 03H/3 is the line flux ratio [Om]5007/H/3 and N2Ha is 
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Figure 17. The difference between An a estimated using equation 12 and 
directly using the Balmer decrement, as a function of a metallicity (from 
Equation 13). The results show that the empirical calibration derived is able 
to use [On]/Ho! line ratios as dust extinction probe without any signifi- 
cant bias for a relatively wide range of metallicities/abundances. A com- 
parison with what would be obtained using the Calzetti law directly is also 
shown, highlighting that the polynomial fit is able to provide a calibration 
of [On]/Ho! which is less affected by metallicities. Solar metallicity corre- 
sponds to 12+log 10 of 8.66; the figure presents metallicities varying from 
w 0.2 to w 2.2 solar metallicities. Note that while Equation 12 provides re- 
liable results (without any significant systematics/biases) for galaxies with 
abundances (12+log 10) from 8.2 to 8.75, it will overestimate the median 
Ah q by about +0.2 mag for very sub-solar and very super-solar metallicity 
galaxies. 



the line flux ratio [Nll]6584/Ha. This indicator has the main ad- 
vantages of i) using emission lines which have very similar wave- 
lengths, thus being essentially independent of dust attenuation and 
ii) having a unique metallicity for each line flux ratio. Figure [T7] 
shows the difference between Ah« computed with the empirical 
[Oll]/Ha and Ah q estimated directly from the Balmer decrement, 
as a function of metallicity. For comparison, the Calzetti law pre- 
diction is also shown, emphasizing that the mismatch between the 
latter law and the observational data is mostly due to the the ef- 
fect of metallicity on the [Oll]/Ha ratio, which is not taken into 
account by Calzetti, but is incorporated in the empirical calibration 
presented in this work. The results suggest that even if galaxies 
at z = 1.47 have different metallicities from those in SDSS, no 
significant offset (within the scatter, ±0.3 mag) is expected when 
estimating Ah« from [Oll]/Ha for a wide range of metallicities. 

For the remaining of the analysis, Ah« is computed as in equa- 
tion 12, both for the z = 1.47 and the SDSS samples. It should 
be noted that the qualitative and quantitative results remain un- 
changed if Balmer decrements are used instead to estimate Ah« 
for the SDSS galaxies, and that qualitative results also remain un- 
changed if the Calzetti law/best linear fit is used instead of the 
polynomial fit. The sample of Ha emitters at z = 0.1 presents 
A H c = 0.97 ± 0.42, while the sample of Ha emitters at z = 1.47 
presents Ah« = 1.0 ± 0.6. 



5.3 



0.1 



[Oll]-Ha luminosity correlation at z = 1.47 and z 

Figure [9] presents the distribution of [Oll]/Ha line ratios for the 
z = 1.47 sample. While it reveals a relatively wide range within 
the sample, (« 0.08 — 1.2), it also shows that down to the Ha 
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Figure 18. A comparison between observed [On] and Hen luminosities (not 
corrected for dust-extinction) for samples at z ~ 0.1 and z ~ 1.5. The 
correlation between observed luminosities seem to evolve very little in the 
last ~ 9 Gyrs, revealing a typical 1 mag of extinction at Ha for the probed 
luminosities. Lines of constant extinction are also shown - note that these 
are not evenly spaced due to the use of the non-linear relation (Equation 
12). 



flux limit, the line ratio distribution peaks at « 0.5 (a bit above 
the median value, « 0.45). In the SDSS-derived sample, the line 
ratios show a similar range from 0.1 — 1 (c.f. [Hopkins et al.| 2003 ) 
and the median observed line ratio is found to peak at w 0.4. Thus, 
even though the Ha and [Oil] luminosities in the z = 1.47 sample 
are much higher than those probed locally, the typical line ratio, if 
anything, is higher, indicating lower extinction. 

This is investigated in more detail in Figure [18] which shows 
how the Ha and [Oil] luminosities (not corrected for extinction) 
correlate over the range of luminosities probed by SDSS at z ~ 0.1 
and HiZELS at z = 1.47. Lines of constant extinction (in Ha) 
are also shown. It is noteworthy that at a given observed Ha lumi- 
nosity, the median [Oil] luminosities are slightly higher (indicating 
lower median extinction) at the higher redshift. Both samples show 
weak trends for increasing extinction with increasing luminosity. 
Furthermore, it should be noted that the simple assumption that 
both samples present a typical constant extinction of 1 magnitude 
is a relatively good approximation. 

Hop kins et al.| < |2001) ) find that there is a correlation between 
SFR and Ana, and so argue that it is possible to estimate (statis- 
tical) dust extinction corrections based on observed SFRs, partic- 
ularly for Ha (but also for the [Oil] -derived SFRs); their relation 
has been used to apply statistical corrections to the observed Ha 
luminosities in many studies at low and high redshift. On the other 
hand, the Hopkins SFR- Ah a relation, whilst derived with a rela- 
tively small sample in the local Universe, seems to be roughly valid 
at z ~ 0.84, as found by Garn et al. (2010 ) (which presents a sim- 
ilar analysis to Hopkins et al., i.e. comparing mid-infrared SFRs 
with Ha), although it slightly overestimates the amount of dust ex- 
tinction of the sample. Nevertheless, it is unclear whether a similar 
result can be found when studying Ana as a function of SFRs at 
z — 1.47, or if the result from [Hopkins et aL](2001) can be re- 
covered when using emission-line ratios (e.g. Balmer decrement, 
or the [Oll]/Ha calibration) at z ~ 0.1. 

Figure [l9]presents the results of investigating the dependence 
of dust-extinction on observed Ha SFR for both samples. The re- 
sults show a relatively weak correlation between Ah a and observed 



' Ha: + [Oil] (2=0.1) ■ 

^7 , Median (z=0.1) i 

. | Ha:^[0 + II] + (Z=1.47) | 

Median (z=L47) ■ 



1 1 1 1 1 Hopkins relation r ■ u 




12 5 10 25 50 

Observed Ha SFR (M,.,yr l ) 



Ha + [Oil] (z=0.1) / 
V Median (z=0.1) / , 

■ Ho + [Oil] (z=1.47) / +v / 
+ Median (z=1.47) /^^V 




4.0 
3.5 
3.0 
"B 2.5 

CO 

-§-2.0 

a 
X 

< 1.5 
1.0 
0.5 



4.0 
3.5 

3.0 
B2.5 

-§-2.0 

a 

< 1.5 
1.0 
0.5 

1 2 5 10 25 50 100 250 750 
Extinction-corrected Ha SFR {M yr _1 ) 

Figure 19. The relation between dust extinction and star formation rates 
(uncorrected and corrected for dust extinction) based on the Ha: emission 
lines for both z = 1.47 and z = 0.1. The top panel shows the relation be- 
tween dust extinction at Hen and observed Hen SFRs. The figure also shows 
the |Hopkins et al.|{2001) relation; this consistently over-predicts the dust 
extinction correction of the z = 1.47 sample by ~ 0.5 mag and, although 
it agrees better with the z = 0.1 sample, it over-predicts the slope of the 
correlation. The bottom panel presents the relation between An a and dust 
extinction-corrected SFRs for both redshifts. Note that because both sam- 
ples are Ho; selected, there are strong selection biases affecting the Hoi anal- 
ysis (dashed lines, resulting from the approximate selection limit in the two 
samples, and from using An a to correct SFRs for extinction), but the trends 
are still recovered even when accounting for such biases. The results show 
that even though the trends are similar at both z = 0.1 and z = 1.47, the 
normalisation is different, and the trends have evolved between z ~ 1.5 
and z ~ 0. 



Ha SFRs at both z = 0.1 and z = 1.47; an offset in the median 
extinction for a given SFR between both epochs is also found. 

The figure also shows the Hopkins et al. ( 2001 ) relation be- 
tween Ah« and observed SFRs, with the results showing that 
it consistently over-predicts the dust extinction correction of the 
z — 1.47 sample by ~ 0.5 mag (~ 1.6 x) and, although it agrees 
better with the z — 0.1 sample (in the normalisation), it over- 
predicts the slope of the correlation. Note that Ah q at z — 1.47 
are potentially overestimated for extremely metal-poor or metal- 
rich galaxies, so the offset is robust and, if anything, is under- 
estimated. The SDSS results remain completely unchanged when 
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using the Balmer decrement Ah«. Therefore the use of the Hop- 
kins et al. relation to correct observed Ha SFRs as a function 
of observed SFRs/luminosity results in a clear overestimation of 
the dust-extinction correction at z = 1.47; re-normalising it by 
^0.5 mag is able to solve this. 

Nevertheless, there seems to be a correlation between SFR and 
Ahc* at z ~ 1.5, which is more clearly revealed after correcting 
SFRs for dust extinction, as can be seen in Figure [19] It should 
however be noted that such relation is in part a result of a bias (as 
indicated in Figure [19]), as for a given flux limit and a (wide) distri- 
bution of dust extinction corrections, one easily recovers a relation 
between extinction and corrected SFRs (as correcting SFR for ex- 
tinction relies on Ana)- The relation between extinction and SFR 
is clear at the highest Ha luminosities where such selection biases 
are negligibly small. The offset between the Ah« vs corrected-SFR 
relations between z = 0.1 and z — 1.47 is also clear at these high 
SFRs, and is still recovered even when a common Ha luminosity 
limit is applied to both samples. These results indicate that although 
there does appear to be a relationship between dust extinction and 
SFR, this relation appears to evolve with redshift and should not 
be used as a reliable way of estimating statistical dust extinction 
corrections for samples of galaxies at different redshifts. 

5.4 Mass as a dust-extinction indicator 



Recently, |Garn & Best (2010 ) performed a detailed investigation of 
the correlations between dust extinction (Ah«) and several galaxy 
properties (e.g. metallicity, star formation rate, stellar mass) using 
large SDSS samples. The authors find that although Ah q roughly 
correlates with many galaxy properties, stellar mass seems to be 
the main predictor of dust extinction in the local Universe (see also 
Gilb ank et al.||2010[ for a similar analysis). The authors derive a 
polynomial fit to the observed trend, which can be used to estimate 
dust extinction corrections for galaxies with a given stellar mass in 
the local Universe. Nonetheless, so far no study has been conducted 
in order to test whether such relation exists at high redshift and 
whether it evolves significantly. 

5.4.1 Estimating stellar masses at z — 1.47 

In order to investigate any potential correlation between dust ex- 
tinction and stellar mass at z ~ 1.5, stellar masses are obtained for 
the entire z — 1.47 sample, following the methodology fully de- 
scribed in Sobral et al. ( 201 1) >. Very briefly, the multi- wavelength 
data available for the z = 1.47 sources are used to perform a full 
SED x 2 fit with a range of models - normalised to one solar mass 
- to each galaxy; the stellar-mass is the factor needed to re-scale 
the luminosities in all bands from the best model to match the ob- 
served data. As in Sobral et al. ( 201 1), the SED templates are gen- 
erated with the stellar population synthesis package developed by 
Bruzual & Chariot (2003), but the models are drawn from B ruzual| 
(2007]. SEDs are produced assuming a universal initial mass func- 
tion (IMF) from Chabrier ( 2003) and an exponentially declining 
star formation history with the form e _t//r , with r in the range 
0.1 Gyrs to 10 Gyrs. The SEDs were generated for a logarithmic 
grid of 220 ages (from 0.1 Myr to 4.3 Gyr - the maximum age at 
z = 1.47). Dust extinction was applied to the templates using the 
Calzetti et aT]([2000) law with E(B - V) in the range to 0.8 (in 
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Figure 20. The relation between dust extinction (Ah q ) and stellar mass 
for the samples at z = 0.1 and z = 1.47. Individual galaxies from both 
samples and the median are presented. The relation derived by Garn & Best 
( 2010) is also shown, which is found to describe the relation very well for 
z = 0.0 — 1.5 at least for low and moderate masses, indicating that stellar 
mass is likely a fundamental property for predicting dust-extinction cor- 
rections. The most massive star-forming galaxies (log Mass> 10.75) may 
have different extinction properties at z = 0.1 and z = 1.5, although that 
may well be a result of the small number of such galaxies at z = 1.5. 



direct comparison with the observed total fluxes. For all except 5 
sources (1 with no U band data, and 4 not detected in any IRAC 
band), 16 bands are used, spanning from the CFHT U band in the 
near-ultra-violet to the 4 IRAC bands (M. Cirasuolo et al., in prep.). 
The appropriate corrections (c.f. Sobral et al. 2011) are applied to 
obtain total fluxes in each band. 

Stellar mass estimates of each individual source are found to 
be affected by a la error (from the multi-dimension x 2 distribu- 
tion) of ~ 0.30 dex, which results from degeneracies between the 
star formation time-scale r, age, extinction and, to a smaller extent, 
metallicity. As the analysis uses a Chabrier IMF, stellar masses are 
directly comparable with SDSS masses (which are the ones used 
by |Garn & BestpOlO) , but note that there is a systematic offset 
when compared to Salpeter (see |Sobral et al.| 2011). It should also 
be noted that the E(B — V) from the best fits correlate well with 
the Ah a determined for each individual galaxy. 



5.4.2 Mass-Extinction relation 

The sample of z — 1.47 emitters presents a median stellar mass 
of 1O 9 9 M . Figure 



20 



steps of 0.1). The models are generated with a logarithmic grid of 
6 different metallicities, from sub-solar to super-solar metallicity. 
It is assumed that all Ha emitters are at z = 1.47 and the com- 
plete filter profiles are convolved with the generated SEDs for a 



shows the observed relation between Ah a 
and stellar mass, for both z = 0.1 and z — 1.47 samples, together 
with the |Garn '& Best ( 2010) relation. The results reveal that not 
only is there a correlation between stellar mass and dust extinction 
at z = 1.47, just like the one at z ~ 0, but, even more importantly, 
that the Garn & Best ('2010) relation seems to be valid at least up to 
z ~ 1.5 as a dust extinction estimator for most masses. As shown 
above, this is in contrast with SFR-dependent extinction correc- 
tions, which must at least be re-normalised when being applied to 
z ~ 1.5 or 2 ~ 0, and provides an important insight into what 
is important in determining the dust properties of galaxies. Table [4] 
presents the Garn & Best ( 2010) relation for predicting Ah« as a 
function of stellar mass, which can be applied, at least within the 
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Figure 21. The relation between dust extinction and rest-frame u — z colour 
for the samples at z = 0.1 (SDSS) and z = 1.47. The data are used to 
derive an empirical relation which is valid for Hoi emitters both at z = 0.1 
and z = 1.47 and predicts median dust extinction corrections (An a ) based 
onu — z colours. Table [4] presents the best-fit polynomial relations based 
on other rest-frame colours which also correlate well with An a . 



studied range of masses, to derive statistical corrections for samples 
of galaxies up to z ~ 1.5. 

Although the relation between stellar mass and Ah« seems to 
hold across redshift for most masses probed, there is tentative evi- 
dence of an offset at the highest masses, in the sense that the most 
massive z = 1.47 star-forming galaxies appear to be affected by 
significantly less dust extinction that those with comparable masses 
at z = 0.1. The sample at high masses remains relatively small, 
however, and also since these massive star-forming galaxies are all 
of very low Ha EWs (see Sobral et al. 2011), there is the possi- 
bility that selection effects are driving this result. Until this can be 
investigated with an improved sample, however, caution should be 
taken in applying the mass-extinction relation at high redshifts for 
masses above ~ io 10 75 M . 



Table 4. Predicting median dust-extinction (Ah q , mag) corrections using 
rest-frame colours and mass for galaxies at z ~ — 1.5. Best- fit relations 
are based on 3rd order polynomials: y = Ax 3 + Bx 2 + Cx + D. Lower 
and upper limits of validity of the fits are also presented. Note that the pho- 
tometry is in the AB system and using SDSS filters and that Mass is given 
in units of log(M0 /10 10 ). The scatter of the data relative to the best fit is 
also given. 



Property 


A 


B 


C 


D 


Validity 


Scatter 


Mass 


-0.09 


0.11 


0.77 


0.91 


[-0.4,0.3] 


0.33 


(u - g) 


-1.31 


4.59 


-4.15 


1.68 


[0.4,1.5] 


0.35 


(g - r) 


-2.63 


5.00 


-0.78 


0.51 


[0.0,0.8] 


0.28 


(r - i) 


-27.29 


29.83 


-7.64 


1.10 


[0.1,0.55] 


0.35 


(u - z) 


-0.092 


0.671 


-0.952 


0.875 


[0.5,3.2] 


0.30 



profiles broadly match u and z SDSS filters and a simple statistical 
correction of +0.1 j^]on the (z — K) colour (which accounts for a 
small K correction and the differences in the filter profiles) is able 
to match the colours very well. 

The results presented in Figure |2l] show that there is a sig- 
nificant correlation between galaxy colour and dust-extinction and 
suggest that, despite galaxies in the sample at z — 1.47 being bluer 
(on average), a single relation seems to hold across epochs (at least 
out to 2 ~ 1.5). Indeed, a simple polynomial fit to the median 
extinction for galaxies with a given rest-frame u — z colour (« ob- 
served (z' — K) + 0.15 colour at z ~ 1.5) is valid at both z — 0.1 
and z = 1.47, and is given by: 

Arc = -0.092(u-z) 3 +0.671(u-z) 2 -0.952(u-z)+0.875(l4) 

Relations between Ah a and various other optical rest-frame 
colours are also investigated; these can be a valuable tool to es- 
timate dust-extinction of galaxy populations at different redshfits 
where only a simple colour is available. Table [4] presents the best 
fits to the data that are valid at least up to z ~ 1.5, together with 
the limits within the relations are valid. The scatter is also quanti- 
fied for each fit (see Table [4), revealing that the relations provide 
good fits of comparable quality to the mass-extinction relation. 



5.5 Predicting A Hc * dust-extinction with colours 

Stellar mass seems to be a relatively reliable (and important) non- 
evolving dust-extinction predictor which can therefore be applied 
both in the local Universe and at higher redshift (z ~ 1.5), at 
least within a range a masses. However, computing reliable stellar 
masses requires obtaining a significant number of observations at 
different wavelengths, relies on the validity of the stellar population 
models used to compute mass to light ratios, and it is, in general, a 
difficult quantity to estimate, particularly at high redshift. There is 
therefore the need to investigate more direct observables as a way 
to predict the median dust extinction of a population of galaxies 
that could be applied out to z ~ 1.5. Rest-frame optical colours 
are expected to correlate with dust extinction, and are therefore in- 
vestigated and calibrated as dust-extinction tracers. All colours pre- 
sented in this section are given in the AB system and are rest-frame 
colours computed using the SDSS filters. 

Figure [21] presents Ah« as a function of rest-frame (u — z) 
for both z — 0.1 and z — 1.47. For the z — 1.47, the observed 
{z - K) colour^ can be used, as the redshifted z and K filter 

13 2' -band data from Suprime-cam/Subaru; K from WFCAM/UKIRT. 



5.6 Discussion of emission-line ratios 

The mean dust extinction properties of the sample of moderately 
star-forming galaxies at z ~ 1.5 seem to be very similar to those 
in the local Universe (a simple 1 mag of extinction at Ha for the 
entire population of z ~ 0.1 and z ~ 1.5 galaxies is a relatively 
good approximation, but with a scatter of 0.3 mag), as a whole, 
even though modest-SFR galaxies at z ~ 1.5 seem to be slightly 
less extinguished. Even more interesting is the fact that dust ex- 
tinction presents the same dependence on stellar mass in the last 
9 Gyrs, at least for star-forming galaxies with low and moderate 
stellar masses. In contrast, while dust extinction correlates with 
SFRs at both z — and z ~ 1.5, the normalization of the rela- 
tions clearly evolves, with differences of ^ 0.5 mag in Ha for the 
same (corrected) SFR. 

As extinction-corrected SFRs correlate reasonably well with 



14 This statistical correction was computed by using the range of best-fit 
SED models to the z=1.47 sources and measuring the (z f — K) colours 
with the z' and K filters on Subaru and UKIRT, respectively, following by 
measuring the (u — z) rest-frame colours of the same sources, using the u 
and z SDSS filters. 
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stellar mass and dust extinction also correlates with (corrected) 
SFRs, it is possible that the relation simply evolves as sSFRs 
evolve. Physically, the normalization of the relation could be driven 
by the gas reservoirs in galaxies; allowing them to reach much 
higher SFRs at z ~ 1.5 than locally, for a fixed stellar mass. This 
conclusion is in line with Garn & Best (2010) and has important 
consequences towards understanding galaxy evolution in the last 9 
Gyrs and how little dust properties seem to have changed in galax- 
ies with modest SFRs. 



6 CONCLUSIONS 

This paper presented the results from the first panoramic matched 
Ha+[Oll] dual narrow-band survey at z ~ 1.5. This is a very effec- 
tive way of assembling large and robust samples of Ha and [Oil] 
emitters at z = 1.47. It provides a large, robust sample of Ha emit- 
ters at z = 1.47, together with a large sample (~ 1400) of [Oil] 
emitters at the same redshift. The survey has allowed for the first 
statistical direct comparison of Ha and [Oil] emitters at z ~ 1.5 
and a direct comparison with an equivalent sample in the local Uni- 
verse to look for evolution. The main results are: 

• The well-defined samples of emitters were used to compute 
the Ha and [Oil] luminosity functions at the same redshift. For the 
Ha luminosity function at z — 1.47, the best-fit Schechter func- 
tion parameters are: logL^^ = 42.5 =b 0.2 ergs -1 , log^c = 
-2.4 ± 0.3 Mpc -3 and a H « = -1.6 ± 0.4, while for the [Oil] 
luminosity function at the same redshift the best-fit parameters 
are: logLf OII] = 41.71 ± 0.09 ergs -1 , log0f OII] = -2.01 ± 
0.10 Mpc -3 and a [ou] = -0.9 ± 0.2. 

• Both Ha and [Oil] luminosity functions show a strong and 
consistent evolution in 0* and L* from z ~ to z ~ 1 and a 
continued L* evolution to z ~ 1.5 and beyond. By combining 
the results with other HiZELS measurements and other estimates 
from the literature, our understanding of the star-formation his- 
tory of the Universe is improved. Using a single well-calibrated 
indicator, the star-formation rate density is shown to rapidly in- 
crease out to z ~ 1, and to probably continue rising (although 
much more weakly) out to z ~ 2, due to the steep faint-end 
now measured at z = 2.2. At z = 1.47, the Ha analysis yields 
p SFR = 0.16 ± 0.05 M yr -1 Mpc -3 , while the [Oil] analysis 
yields /?sfr = 0.17 ± 0.04 M yr -1 Mpc -3 , in excellent agree- 
ment. 

• By using SDSS, the [Oll]/Ha line fraction is calibrated as a 
dust-extinction probe against the B aimer decrement. The relation 
is shown to be accurate within ^0.3 dex. 

• Ha and [Oil] luminosities correlate well at z = 1.47, simi- 
larly to those at z ~ 0.1, but the sources at higher redshifts appear 
less dust extinguished for a given observed Ha luminosity. A rela- 
tively weak correlation between observed SFR and dust extinction 
is found for both z = 0.1 and z = 1.47, but with a different nor- 
malisation. It is also shown that the Hopkins relation consistently 
over-predicts dust extinction corrections for z ~ 1.5 by ~ 0.5 mag 
in Ha. 

• Stellar mass is shown to be a good dust-extinction predictor, 
at least for low and moderate mass galaxies, with the relation be- 
tween dust extinction and mass being the same in the last 9 Gyrs 
for such star-forming galaxies. The relation between mass and dust 
extinction from Garn & Best (2010) is shown to be fully valid with 
no evolution at z = 1.47. 

• Optical or UV colours are shown to be a simple observable 



extinction predictor which can be applied for z ~ — 1.5 star- 
forming galaxies; the best-fit relations based on several colours are 
derived and presented. 

The results presented in this paper contribute to our under- 
standing of the nature and evolution of star-forming galaxies at the 
likely peak of the star-formation history of the Universe, particu- 
larly by showing an effective and clean way of selecting large sam- 
ples of these galaxies and by investigating, for the first time, the 
evolution using both the Ha and [Oil] emission lines at z ~ 1.5 
at the same time. The results show a very good agreement between 
the Ha and [Oil] view. Moreover the results reveal that the typi- 
cal necessary extinction corrections for the probed Ha luminosi- 
ties at z = 1.5 is Aua ~ 1.0 (but with a scatter of « 0.3 mag); 
this is what has been found for a range of luminosities in the local 
Universe, therefore revealing no significant evolution for moderate 
SFR galaxies (if anything, sources at higher redshift are less ex- 
tinguished). Extinction corrections blindly applied as a function of 
Ha luminosity (e.g. using the local Hopkins relation) over-predict 
the dust extinction correction for the z — 1.47 sample and would 
lead to significant biases. 
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APPENDIX A: COMPLETENESS CORRECTIONS 

Section |4~2| outlines the procedure which has been followed in or- 
der to address such issues and further details are given in S09. Here 
details are given on the effect of the changes in methodology since 
S09, by addressing how completeness estimations change with dif- 
ferent input samples for the colour- selection analysis. 

In S09, the entire range of sources detected in the narrow-band 
imaging was used to add line fluxes and then study the recovery 
rate. This can, nonetheless, be improved (as detailed in Section 
|4.2| ), to compute completeness corrections specifically for z ~ 1.5 
sources. Here, that is done by i) excluding all stars (since they are 
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Figure Al. Observed H band magnitude (R rest-frame at z = 1.47) dis- 
tributions for the two samples used in the completeness simulations. Note 
the difference between the samples at bright magnitudes: while the sample 
of likely z > 1 galaxies presents a strong peak at faint magnitudes and a 
relatively sharp decline for brighter magnitudes, the sample containing all 
detections presents a much broader distribution, including a much shallower 
decay at bright magnitudes. 
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Figure A2. A study of the completeness fraction (defined as the fraction 
of sources with a given flux recovered by the selection against the actual 
number of source with that flux) using two different input samples (all de- 
tections which are non-emitters and only those consistent with being z > 1 
sources). 



not real z ~ 1.5 galaxies and they do not have their properties) and 
ii) excluding sources with colours which clearly place them clearly 
at z < 1. By taking this approach, the observed H band magnitude 
distribution is quite different from that of all detections, as can be 
seen in Figure [AT] Indeed, by using the entire population of detec- 
tions, the number of bright sources with low equivalent widths is 
overestimated, and thus an overall lower completeness fraction is 
derived. Figure |A2] compares the completeness fraction as a func- 
tion of input line flux for simulations using all sources and for those 
using only z > 1 galaxies (rejecting potential stars and using the 
BRiK colour-colour selection). 



